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The accuracy of spacecraft attitude control using magnetic actuators only is low and on the order of 0.4-5 
degrees. The key reason is that the magnetic torque is two-dimensional and it is only in the plane 
perpendicular to the magnetic field vector. In this paper novel attitude control algorithms using the 
combination of magnetic actuators with Reaction Wheel Assembles (RWAs) or other types of actuators, such 
as thrusters, are presented. The combination of magnetic actuators with one or two RWAs aligned with 
different body axis expands the two-dimensional control torque to three-dimensional. The algorithms can 
guarantee the spacecraft attitude and rates to track the commanded attitude precisely. A design example is 
presented for Nadir pointing, pitch and yaw maneuvers. The results show that precise attitude tracking can 
be reached and the attitude control accuracy is comparable with RWAs based attitude control. The 
algorithms are also useful for the RWAs based attitude control. When there are only one or two workable 
RWAs due to RWA failures, the attitude control system can switch to the control algorithms for the 
combined magnetic actuators with the RWAs without going to the safe mode and the control accuracy can be 
maintained. 


Nomenclature 


bit) 

= Earth’s magnetic field vect 

D{t),k{t) 

= control gain matrices 

J 

= inertia matrix 

m 

= vector of magnetic dipoles 

Sq 

= attitude error quaternion 

q 

= quaternion 

T 

= control torque 

(O 

= angular velocity 

Q(<y) 

= skew- symmetric matrix 


I. Introduction 

Magnetic torque bars as attitude control actuators have been used and studied for many years. The typical uses of 
magnetic actuators in attitude control are for Reaction Wheel Assembles (RWAs) de-saturation, and angular- 
momentum and nutation control of spinning and momentum-biased spacecraft. There are some studies Ref. 1-6 to 
use the magnetic torque bars only to control the spacecraft attitude. Global stable attitude control laws have been 
investigated for both inertial and earth pointing in Ref. 1-2. Periodic control algorithms have been developed using 
an asymptotic linear quadratic regulator technique for a Nadir-pointing Ref.3. A pseudospectral control law has been 
applied to the magnetic attitude control of satellites in elliptic low earth orbits in Ref. 4. 

The main drawback of the magnetic actuator only attitude control is that the control accuracy is low and on the 
order of 0.4-5 degrees. The key reason is that the magnetic torque is two-dimensional and it is only in the plane 
perpendicular to the magnetic field vector. The spacecraft cannot be controlled precisely in three-dimensional space 
using only magnetic actuators. In order to improve the control accuracy, control algorithms for a combination of 
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magnetic actuators with RWAs or other types of actuators are presented in this paper. The combination of magnetic 
actuators with one or two RWAs aligned with different body axis expands the two-dimensional control torque to 
three-dimensional. The attitude control accuracy for the combined actuators using the control algorithms proposed in 
this paper is comparable with RWAs based attitude control. 

The control algorithms of the combined actuators are also useful for the RWAs based precise attitude control. 
When there are only two or one workable RWAs due to RWA failures, the attitude control system can switch to the 
control algorithms for the combined magnetic actuators with the RWAs presented in this paper without going to the 
safe mode. The control accuracy can be maintained. 

This paper is organized as follows. Section II presents the spacecraft dynamics, the kinematics equations, and the 
control matrix produced by the magnetic actuators. The controllability of magnetic actuators only is investigated. 
The necessity of the combination of magnetic actuators with RWAs to control the spacecraft attitude precisely is 
revealed. The combined control matrices with one RWA along different body axis are derived in Section III. The 
attitude control of magnetic actuators with two RWAs along different body axis is discussed in Section IV. Attitude 
tracking algorithms are described in Section V and a design example is given in Section VI to demonstrate the 
performance of the control system. 


II. Problem Statement 


The spacecraft dynamics in the spacecraft body frame can be expressed as 

J6) = —co x (Jco) + T mtb +T rwa + T dist 
= Cl(co)Jco + T mtb +T rwa + T dist 


( 1 ) 


where J is the inertia matrix. 0 ) = (co x , 


CO,. 


CO. 


M 5 is the angular velocity with respect to (wrt) an inertial 


frame expressed in the body frame. T mtb and T RWA are the MTB and RWA control torques. T dist is the 

external disturbances torque due to gravity gradient, aerodynamic force, and solar pressure. Cl(co) is a skew- 
symmetric matrix defined by 


Q(co ) 


0 

-co. 


— co„ 




CO. 

0 co x 

-co x 0 

Using quaternions, the kinematic equation of motion can be described by 


V 


( 2 ) 


J 


^Q(co)q + ^q 4 co 


< 1 4 


1 


( 3 ) 


T 

co q 


The spacecraft attitude is determined by its quaternion q , which represents the attitude of the spacecraft wrt the 
inertial frame. Both q and q 4 of the quaternion q are defined as 


( 4 ) 





r e x sin(^/ 2) ' 

q = 

q 2 

= 

e 2 sin(^/2) 


V^3 J 


v e 3 sin(^/2) y 


q 4 =cos(^/2) 
q=\ 


r q ' 


W4 


( 5 ) 

(6) 


In Equation (4), e l , e 2 , and e 3 are the components of the rotation axis unit vector expressed in the inertial 
frame; (j) is the rotation angle. 
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A set of three magnetic coils aligned with the body frame axes generates the magnetic attitude control torque, 
which can be written as 

T mtb =rnxb(t) = nip ( t)}n (7) 

where m- (m x , m v , m. is the vector of magnetic dipoles for the three coils. b(t) e R 3 is the Earth’s 

magnetic field vector expressed in the body frame, which is changing all the time and depends on the spacecraft 
position and attitude. Note that the magnitude of b (t) is not zero. The following unit vector in the same direction of 
b(t) is introduced as 


b(t) = P^ = (b x (t), b v (t), bXtjf- 

p(0| 

Obviously the unit vector bit) satisfies the following equation 

b 2 x (t)+b 2 y (t)+b 2 ,(t) = 1. 

Consider following form of the vector of magnetic dipoles 


m 


1 

b(t) 


-n T (b(t)\i 


( 8 ) 

(9) 

( 10 ) 


where U = [u x , U y , U^f e R } is a new control vector. The magnetic control torque T mtb for the magnetic 


dipoles is 


where 


t mtb =n{b(t))n T (b(t))u 
= r(Zt(0)w 


(ii) 


r(K0)=^(0M6(0) 


f ul 


b;(t) + b:(t) -b x (t)b y {t) -bJj)bJt) 
-b x (t)b v (t) b;{t)+b:{t) -b v (t)bAt) 


( 12 ) 


-b x (t)b z (t) -b y {t)b z {t) bl(t) + bl(t) J 

We call r(6(?)) the magnetic control matrix. The control capability and spacecraft dynamic behavior using 
magnetic actuators are determined by r(b(t)). The following are some properties of V(b{t)) . 

• r(ft(?))is symmetric and semi-positive definite matrix. It is easy to show the determinant of r(&(?)) is 
zero. 

det|r(*(0) = 0 (13) 

• r(/>(?)) always has one zero eigenvalue and two positive real eigenvalues. The eigenvector corresponding 
to the zero eigenvalue is b(t) given by 

r(6(o)KO = o. (i4) 


• If T MTBj t is the desired magnetic control torque, there is a control vector U only if T mtb ^ ( is 

perpendicular to bit ) . Therefore, in general, if T mtf is not perpendicular to bit) , there is no U to 
realize the desired magnetic control torque. 

The properties of the magnetic control matrix show that the magnetic control torque T mtb is always in the plane, 
which is perpendicular to bit) . The magnetic control torque is only two-dimensional. In order to expand the control 
to three-dimensional, additional actuators need to be added, which can be thrusters or RWAs. In following RWAs 
are considered, but the analysis and principals are applicable to the thrusters. 
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III. Combination Of Magnetic Actuators With One RWA 

If one RWA is employed, there are three configurations, in which the wheel can be aligned with body frame X, 
or y , or z axis. 


Configuration#!: RWA along body X axis 

The control torque T RWA contributed by the RWA along body X axis is given by 

T rwa = («*, 0, Of . 

The combined control torque from the magnetic actuator and the RWA can be written as 

T MTB +T RWA= r { b (t))u + {u x , 0, Of 
= A A b (t))u 

where 



f 

' 1 0 

o v 

\ 

A,(6(/))= 

r(fe( 0 )+ 

0 0 

0 

u 


V 

v 0 0 

Ox 

/ 


r \ +b 2 y {t)+b 2 z {t) 

~b x {t)b y (t) -b x (t)b z «A 

= 

-b x (t)b y (t) 

b 2 x (t) + b 2 (t) - /> (/)/> (/) 


K -b x (tyb,(t) 

~b y {t)b z {t) b 2 (t ) + b 2 (t f 


A x (b(t)) is the combined control matrix with RWA along body X axis. 


(15) 

(16) 


(17) 


Configuration #2: RWA along body y axis 

For the RWA along body y axis the combined control torque is given by 

T MTB +T RWA= T {b(t))u + (0, U y , Of 

= \{b(t))u 

where 


A>(0) 


'b 2 y {t) + b]{t) 
-b x (t)b v (t) 

K ~b x m(t) 


~b x (t)b y (t) 
l + b 2 (t) + b 2 (t) 
~b y {t)b z {t) 


-b x (t)b z {tA 

~b y (t)b z (t) 
b x (t) + b;(t)y 


(18) 


(19) 


Configuration #3: RWA along body z axis 

Similar to other two configurations, the combined control torque can be written as 


T +T 

1 MTB ' 1 RWA 


A >(0) = 


= r(6(0f + (o, o, 

= A z {b{tj)u 

r b 2 y (t)+b 2 (t) -b x (t)b y (t) 
-b x (t)b v (t) b;(t) + b:(t) 


-b x (t)b z (t) 
b. (/)/> (/) 

\ + b 2 x (t) + b](t) 


( 20 ) 


( 21 ) 


/>,(/)/>(/) ~b y {t)b z {t) 

The combined control matrixes A ; .(Z>(t)) (/ = X, or V, or z) representing the RWA aligned with body 

frame X , or y , or Z axis have the following properties. 

• The combined control matrixes are symmetric and the determinants are given by 

det |A ; . (Z?(r))| = bf(t) (/' = x, or y, or z) . 
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• if b t (t) ^ 0 (/ = x, or y, or z) , A,, (b(t)) (i = x, or V, or z) always has three positive real 

eigenvalues and an inverse. The control torque provided by A .{p(f))u (/ = X, or y, or z) is three- 

dimensional and the spacecraft is fully controlled. 

• If b i (t) = 0 (i - x, or y, or z) , A , (^(0) 0 = *, or jp, or z) has two positive real eigenvalues and 

one zero eigenvalue. The eigenvector corresponding to the zero eigenvalue is b(t ) . At the same time b{t ) 
is perpendicular to the wheel direction. 

For the combination of the magnetic actuators with one wheel the combined control matrixes 
A, to) (i = x, or V, or z) is three-dimensional control except b t (t) = 0 (i = X, or y, or z) , in which it is 
two-dimensional control. 

If the wheel is not aligned with any body frame axis, the combined control matrix should have the same 
properties because the wheel control torque can be transformed to body frame axis by a direction cosine matrix. 


IV. Combination Of Magnetic Actuators With Two RWAs 

If two RWAs are combined with the magnetic actuators, there are three configurations of two wheels aligned 
with the body frame axis as follows. 


Configuration #4: RWAs along body X and y axes 
The control torque T RWA contributed by the RWAs is given by 

TrwA = [ U x’ U y’ o) 7 ^ 

The combined control torque from the magnetic actuator and the RWAs can be written as 

TmTB + r ^RWA = ^{b(t))u + {u x , u y , of 
= A : (b(t))u 

where 

r \ 0 0 ^ 


A>(0) = 


r(6(0)- 


0 1 0 
0 0 0 


r l+b 2 (t)+bl(t) 


'JJ 

~b x {t)b Y (t) 


bAt)h(t) 


-b x (t)b y (t) l+b;(t)+b;(t) /t (/)/» (/) 
~b x (t)b z (t) b(t)b.(t) b 2 x (t) + b 2 (t) 


Configuration #5: RWAs along body V and z axes 

For the RWAs along body y and z axes the combined control torque is given by 

Tmtb + t rwa = r (b(t))ii + (o, u y , u : j 7 
= K(b(t)}i 

where 


f L 2 


A>(0) = 


b 2 y {t) + b 2 z {t) -b x {t)b y (t) 
~b x {t)b(t) l + b 2 (t) + b 2 (t ) 


-b x {t)b z {t) 

~b y (t)b z (t) 


-b x {t)b z {t) ~b(t)b z (t) l + b 2 (t) + b 2 (t ) 


(23) 

(24) 


(25) 


(26) 


(27) 
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Configuration #6: RWAs along body X and z axes 

Similarly, for the RWAs along body X and z axes the combined control torque can be written as 

T MTB + T RWA= r ( h ( t ))lt + {u x , 0 , ItJ 

= A y {b(t))u 

where 


A, (6(0) 


'l + b 2 y (t) + b*(t) -b x (t)b y (t) -b x (t)b z (t) " 

-b x {t)b y {t) b 2 x (t) + b;(t ) -b, (z )b (/) 

v -b x (t)b : (t) -b y (t)b__{t) \ + b 2 x (t) + b 2 y {t)' 


( 28 ) 


(29) 


The combined control matrixes A ; (Zz(Y)) (z = X, or y, or z) have the following properties. 

• The combined control matrixes are symmetric and the determinants are given by 

det | A,. (b(t)) | = 2(l - bf ( t )) (z = x, or y, or z) ■ (30) 

• If |6 i (t)|^l (z = x, or y, or z) , A ,(b(t)) (i = x, or y, or z) always has three positive real 

eigenvalues and an inverse. The control torque provided by A ; (6(f))zz (z = X, or y, or z) is three- 
dimensional and the spacecraft is fully controlled. 

• If |^(0| = 1 (z = x, or y, or z) , A ; .(Zz(t)) (z = x, or y, or z) has two positive real eigenvalues and 

one zero eigenvalue. The eigenvector corresponding to the zero eigenvalue is b(t) . At the same time h(l ) 
is perpendicular to both of the two wheel directions. The combined control matrixes are two-dimensional 
control. 


600 Km Circular Orbit, 50 Deg Inclination, Nadir Pointing 600 Km Circular Orbit, 90 Deg Inclination, Nadir Pointing 




Figure 1 b(t ) with 600km circular orbit and 50 
degrees inclination 


Figure 2 b(t ) with 600km circular orbit and 90 
degrees inclination 
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Unlike the combination of the magnetic actuators with one wheel, the combined control matrixes with two 
wheels are deteriorated to two-dimensional control only if b{t ) is in the direction, which is perpendicular to both of 
the two wheel directions. For two wheels the chance for two-dimensional control is reduced to b(t ) along one 
direction in three-dimensional space from b{t ) in a plane, but additional wheel is required as the costs. In practice 
b(t) is always changing the direction along the spacecraft orbit. There are some orbits and attitudes, in which b(t) 
does not cross the plane of Zr(7) = 0 (/ = X, or y, or z) for one wheel combination configurations. Such an 

orbit is shown in Figure 1, which is a 600km circular orbit with 50 degrees inclination. The spacecraft attitude is 
Nadir pointing. The calculation of b(t ) is based on the dipole model in Appendix FI of the book Wertz (1978). 

As shown in Figure 1, b {t) is not equal to zero for the Nadir pointing spacecraft on the orbit. If the wheel i 


is 


aligned with body V axis, A v {b(t)) always has an inverse. If the wheel is aligned with body X or z axes, b x (t) 
or b z {t) will cross zero periodically. A v (Zi(?)) or A_(6(/)) has no inverse when b x (t) or b,{t ) crosses zero. 

Figure 2 shows b{t ) for a 600km circular orbit with 90 degrees inclination and Nadir pointing. All components 
of b(t) cross zero periodically. b v {t) crosses zero at much slower rate (about every 8 orbits). Flowever, for two 

wheels configurations, if the two wheels are aligned with X and Z axes, A v {b{t)) always has an inverse for both 
orbits because ||& y (Y)|| is always smaller than one. 


V. Attitude Tracking Algorithms and Performance 

For the combination of magnetic actuators with one RWA the spacecraft dynamics (1) becomes 
Jco = H(oj)Joj + A, ( b(t))u (i = x, or y, or z) 


(31) 


Suppose (O c and q c are the commanded attitude rate and the commanded vector part of the attitude quaternion 
to be tracked. Then CO c is the rotation rate of the command frame with respect to (wrt) an inertial frame resolved in 
the command frame. Therefore, CO c and q c satisfy the following kinematic equation of motion 


?c=-| fi (®c)? c +|?4c®c 


Rac=' 


1 T 
2°' q ' 


(32) 


where 


CO, 


=(' 


®jr c ’ ’ 

f 


CO . 


J 




(33) 

(34) 


0 co zc - co y 

Q(co c )= -co, c 0 oj xc 

v ®yc 0 

^c=(<bl c> <hc’ <hJ > 9c=(9c <14 J ( 35 ) 

In order to introduce the control algorithm, the attitude error quaternion §q needs to be defined as follows: 


J 
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^4 c 

<l3c 

-(he 

-he 

'h 

-q 3 c 

<l4c 

(he 

~he 

h 

<hc 

-<hc 

(^ 4 c 

~he 

h 

<h c 

he 

he 

he , 

V^4 


where q c 1 is the inverse of q c , and ( q c 1 • q ) is quaternion multiplication. Note that Sq is still a quaternion, 
which has a physical meaning. If the spacecraft rotates with q c and then rotates Sq further, its attitude is the same 
as if it rotates with q . Therefore, Sq represents the attitude difference between q and q c . Using the attitude error 
quaternion, the control algorithm for the combination of magnetic actuators with one RWA is written as follows: 

u = -A) 1 (b(t)^Q.(co)Joo + D(co -a> c ) + KSq - Jti> c ] (/ = x, or y, or z ) (37) 

where D and K are 3x3 constant gain matrices to be determined. At the orbit point, where 
b i (t) = 0 (/ = X, or V, or z) , A ( . (&(?)) (/ = X, or y, or z) has no inverse, the inverse at the neighbor orbit 
points is used to form the control law in (37). Substituting (37) into (31) gives 

J{d> - cb c ) + D(co - a> c ) + KSq = 0 . ( 38 ) 

Equations ( 38 ) and ( 3 ) form the closed-loop nonlinear dynamic system of the spacecraft attitude with the control 
law ( 37 ). Using dimensionless control gain matrices, Equation ( 38 ) can be rewritten as 

(co- co c ) + D(co- co c ) + KSq = 0 ( 39 ) 

where 

D = J~'D (40) 

K = j'K (41) 

It is not difficult to verify that the following is one of the solutions to Equations ( 38 ) and ( 3 ) 

[(D = CO„ 


[9 = ^ 

However, solution (42) is not an equilibrium point of the system because, in general, CO and q are not constant 
and are the function of time. When (O c = 0 (which implies q c = constant), 0) = (O c = 0 and 
q = q c = constant , solution (42) is an equilibrium point of the system. Solution (42) means that it is possible for 
the spacecraft to achieve the commanded attitude and attitude rate. This control law can be employed for attitude 
tracking under any large and fast maneuver because there is no restriction on (O c and q c . 

When (O c ^ 0 , (42) will not be the solution of the closed-loop system if the term Q (a>)Ja> or Jd) c is not 
included in the control law ( 37 ). The stability analysis will not draw any conclusion concerning the spacecraft’s 
ability to track the commanded attitude. For the control law ( 37 ) without the term J(b c , the solution (42) exists only 

under condition d) c = 0 (co — constant) . This means that the control system can only track the commanded 
attitude when the attitude rate is constant. If (O c = 0 (a rest-to-rest maneuver), the solution (42) exists for the 

control law ( 37 ) without the terms Cl(co)Jco , D(O c , and J(b c . 

If the solution (42) has local asymptotic stability, CO and q (with initial values near the solution) will converge 
to the solution. If the solution (42) has global asymptotic stability, any CO and q will converge to the solution. The 

following theorem shows how to find the constant gain matrices D and K having the globally asymptotic stable 
solution (42). If such a solution exists, the spacecraft can track the commanded attitude with any commanded 
attitude rate. 
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Theorem Suppose that K 1 is symmetric and positive definite, and K l D is positive definite, 
b^t) ^0 (i — X, or y, or z) . Then, the control law (37) ensures that the solution (42) is global asymptotically 

stable and all trajectories CO and q for the closed-loop nonlinear system (38) and (3) converge to CO • and q c as 
time approaches infinity. 

The proof of the theorem can be found in the paper Zhou and Colgren (2005a). The magnetic dipoles vector m 
required by the control law (37) is given by 

1 


m = - 


n 


fit) II 

The corresponding wheel torque is 


(b (/))A-' {b(t)\fl(co)Jco + D(co - co c ) + KSq - Jcb c ] (i = x, or y, or z) ■ 


(43) 


{u x , 0, Of 


(o. 


j 


( wheel along x axis) 
( wheel along y axis) 
( wheel along z axis) 


[( 0 , 0 , 

Control law (37) can be rewritten as 

u = D(t)(co -co c ) + K(t)Sq - a: 1 (b(t))fl(co)Jco + A7 1 ( b(t))jco c 


(44) 


(/ = x, or y, or z) (45) 


where 


D(t) = -A?(b(tj)D (46) 

k(t) = -A-;(b(t))K (47) 

The total control gain matrices D(t ) and K(t ) are time varying. For two RWAs control law (37) becomes 

u = -AT 1 (b(t)jfl(co)J oo + D(co - co c ) + KSq - (/' = x, or y, or z) . (48) 

Equations (38) and (3) form the closed-loop nonlinear dynamic system of the spacecraft attitude with the control 
law (48). Similar to the Theorem, if K ~ 1 is symmetric and positive definite, and K 'D is positive definite. 
||A(?)|| ^1 (7 = x, or y, or z) . Then, the control law (48) ensures that the solution (42) is global stable and all 

trajectories CO and q for the closed-loop nonlinear system (38) and (3) converge to 00 c and ^ c as time approaches 
infinity. The RWAs torques are given by 


f k> °) r 


T = 

1 RWA 


(o, 

0 


y 


(wheel along x, y axes) 
(wheel along y, z axes) ■ 
(wheel along x, z axes) 


m = — 


0, uj 
n require 

Q. T (b (O^ 1 (b(t))[Cl(co)Jco + D(a>-a> c ) + KSq- Jco c ] 


The magnetic dipoles vector ill required by the control law (48) is 

1 

~ |2 

b(t) 


(49) 


(i = x,or y,or z) ■ (50) 


To evaluate the performance of the control laws, 3-2-1 Euler angle rotation, yaw <)'(// . pitch 50 . and roll 5(f) 
corresponding to Sq are employed. Their relations are given by 

dq { = ±(sin(<50 / 2) cos(£0 / 2) cos(£(// 12)- cos (S</> / 2) sin(J0 / 2) sin(J(// / 2)) 

Sq 2 = ±(cos (Scj) / 2) sin(A# / 2) cos(A(// / 2) + sin (S</> / 2) cos^ / 2) sin (Si// 12)) 

Sq } = ±(cos (Scj) / 2) cos(£# / 2) sin( 5y/ 12)- sin( 5cf> / 2) sin(J# / 2) cos(dy/ 12)) 

Sq 4 = ±(cos(S<f> / 2) cos(£# / 2) cos (Si// / 2) + sva(5cj) / 2) sin (SO / 2) sin (Sy/ 12)) 

When the Euler angles are small, they can be approximated by 


( 51 ) 
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The attitude rate error Aft) is given by 


8<j) = 2 c * 7 , 

\ 80 = 2Sq 2 ■ 
8y/ = 2 Sq 3 

A co = co-co r 


( 52 ) 


(53) 


VI. Design Example 

In order to assess the performance of the control law discussed in this paper, a spacecraft attitude control design 
and simulation tool has been developed. The tool has the following functions: 

• Command Generation: generate Nadir pointing, pitch, and yaw maneuver commands 

• Controller: control algorithm proposed in this paper 

• Spacecraft Dynamic Model: spacecraft rigid body dynamics and orbit mechanics model 

• Environmental Models: aero density model (Jacchia Robert), gravity gradient model, geomagnetic model 
(B field), Solar pressure model 

• Solar Ephemeris Model 

The considered spacecraft has following inertia matrix 


J = 


43.8208 

- 7.6332 


101 . 1218 
0.4111 


0.4111 
371 .3937 


(Kg -m y 


It has a solar array with a 1.2758 m surface along body y axis. It flies on a 600 km circular orbit with 50 
degrees inclination. Two configurations of wheel alignments are presented here. The first one is that the wheel is 
aligned with body V axis, in which A (b(t)j always has an inverse for the Nadir pointing spacecraft on the orbit 

as shown in Figure 1. For the second configuration the wheel is aligned with body Z axis, b z (t) will cross zero 
periodically and A_(/t(7)) has no inverse when b z (t) crosses zero. For both configurations the gain matrices D 
and K in the control law are selected as 

3071 1 9466 -D 999A 

^ Kg ■ m 2 


,Kg-nr 



A 2 

.8071 

1 

. 94 66 

-0 

. 3391^ 

D = 

1 

. 94 66 

4 

.4920 

0, 

.0183 


v-° 

.3391 

0 

.0183 

16 

. 4981 y 


f 0. 

5691 

0. 

0865 

-0. 

015A 

K = 

0. 

0865 

0. 

1996 

0. 

0008 


v -o. 

.0151 

0. 

.0008 

0. 

7 331 y 


K 1 and K l D have the following corresponding values 


K ' = 


K l D = 


506.6059 

0 

0 

( 22.5045 
0 
0 


0 

506.6059 

0 

0 


0 

0 

506.6059 


(s 2 ) 


(s)- 


22 . 5045 

0 22.5045^ 

According to the Theorem, the closed-loop system is globally stable, and the system can track the commanded 
attitude. 

Four simulations for the two wheel configurations are presented here. The first simulation results are shown in 
Figure 3 to Figure 5, which is Nadir pointing acquisition for the wheel aligned with body V axis. The spacecraft 
points to the earth initially, but the attitude rate is zero. In order to point to the earth continuously, the pitch rate has 
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to go to the orbit rate. Figure 3 shows that the attitude rate errors converge after about 300 seconds. The maximum 
attitude error during the transient is about 0.9 degrees. After the transient the steady attitude and rate errors are 
comparable with 3-axis wheel control 

5(f) = -7.9 ( arcsec) 

59 = 59 ( arcsec ) 

Sy/ = -0.3 ( arcsec ) 

Aco x = 0.0003 ( arcsec/s ) 

A co y = 0.00001 ( arcsec/s ) 

Acu_=-0.01 ( arcsec/s ) 

The steady attitude errors are caused by the disturbance torques shown in Figure 5. When the disturbance torques 
are not included in the simulation, the steady attitude errors are less than one arcseconds. In order to reduce the 
steady attitude errors, the gain matrices D and K in the control law can be increased, but they are constrained by 
considering additional spacecraft dynamics, such as flexible modes and fuel slosh. The magnetic dipoles usage and 
the wheel torque shown in Figure 4 are reasonable. 

In the second simulation Nadir pointing with wheel along z axis is conducted as shown in Figure 6 and Figure 7. 
The spacecraft points to the earth and rotates with the orbit rate initially. The time period of the simulation is about 
two orbits. There are three times, in which the steady errors have transient. They are corresponding to b,(t) crossing 
zero three times. The absolute value of the minimum eigenvalue of A_(b(t)) is used to measure the singularity. 
When the absolute value is less than 0.0002, the inverse of A_(fo(t)) with the absolute value larger than 0.0002 is 
used to form the control law in the simulation. The singularity causes the transient, but the attitude errors transient is 
small and short. The control torques and the dipoles response significantly, but the maximum usage of the dipoles 
and the wheel torque are reasonable. The attitude control accuracy is comparable with 3-axis wheel control. 

A pitch maneuver starting with Nadir pointing for the wheel aligned with body y axis is conducted in the third 
simulation as shown in Figure 8 to Figure 10. The simulation starts with Nadir pointing attitude and rates. At 500 
seconds the spacecraft begins to pitch up 90 degrees within 20 minutes and hold the 90 degrees pitch angle after the 
maneuver. During the maneuver the attitude and rate errors are small. The maximum pitch angle error is about 300 
arcseconds. The dipoles and wheel torque usages are moderate. The command tracking during the maneuver is 
precise. 

In the last simulation a yaw maneuver for the wheel aligned with body y axis as shown in Figure 1 1 to Figure 
13 is conducted. Different from the pitch maneuver, the Nadir pointing is preserved during the yaw maneuver. The 
yaw maneuver angle is 180 degrees with the duration of 40 minutes. The maneuver starts at 500 seconds. The 

attitude errors around about 2100 seconds have a spike, which is caused by b At ) crossing zero. The maximum 

pitch angle error is about 1.8 degrees, but overall the yaw maneuver tracking is precise. Similar to the second 
simulation, the absolute value of the minimum eigenvalue of A v (b(t)) is used to measure the singularity. When the 

absolute value is less than 0.0002, the inverse of A (b(t)) with the absolute value larger than 0.0002 is used to form 

the control law in the simulation. Further the magnetic dipoles are limited within ± 400 (A —in 2 ) and the wheel 
torque is limited within +0.02 (Nm) to reflect the hardware limitation in real ACS design. The attitude and rate 
errors can be reduced with expansion of the limits for the magnetic dipoles and the wheel control torque. 
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Figure 3 Attitude and rate errors during Nadir pointing acquisition 
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Figure 4 Magnetic dipoles and control torques during Nadir pointing acquisition 
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Figure 5 S/C attitude rates and disturbance torques during Nadir pointing acquisition 
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Figure 6 Attitude and rate errors during Nadir pointing with wheel along z axis 
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Figure 7 Magnetic dipoles and control torques during Nadir pointing with wheel along z axis 
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Figure 8 Attitude and rate errors during pitch maneuver with Nadir pointing 
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Figure 9 Magnetic dipoles and control torques during pitch maneuver with Nadir pointing 
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Figure 10 S/C attitude rates during pitch maneuver with Nadir pointing 


15 

American Institute of Aeronautics and Astronautics 


Yaw Maneuver 180 deg with Nadir Pointing 


Wheel Along y Axis 



Yaw Maneuver 180 deg with Nadir Pointing 
Wheel Along y Axis 







Time (sec) Time (sec) 


Figure 11 Attitude and rate errors during yaw maneuver with Nadir pointing 
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Figure 12 Magnetic dipoles and control torques during yaw maneuver with Nadir pointing 
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Figure 13 S/C attitude rates during yaw maneuver with Nadir pointing 
VII. Conclusions 

Spacecraft attitude control using the combination of magnetic actuators with RWAs or other types of actuators, 
such as thrusters, is presented in this paper. The attitude control algorithms of the combined actuators are derived, 
which can guarantee the spacecraft attitude and rates to track the commanded values precisely. A design example is 
presented for Nadir pointing, 90 degrees pitch and 180 degrees yaw maneuvers. The results show that precise 
attitude tracking can be reached and the attitude control accuracy is comparable with 3 -axis wheel control. 
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